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Abstract 



We examine the evolution of a holographic cosmological model with future event horizon as 
the infrared cut-off and dark matter and dark energy do not evolve independently — there is 
^ ' interaction between them. The basic evolution equations are reduced to an autonomous system 

and corresponding phase space is analyzed. 
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! 1 Introduction 

Recent observational evidences particularly from Type la supernovae and Cosmic Microwave Background(CMB) 
speculate the existence of both gravitating and non-gravitating type of matter. There is a substantial amount 
\ of gravitating matter non-baryonic in nature and is termed as Dark Matter (DM) [1] [51 [31 IH [S] . On the other 

' hand, the non-gravitatng matter, known as Dark Energy(DE) is the mysterious agent for the present phase of 

\ cosmic accelerated expansion. It is only known for certain that DE has huge negative pressure(comparable to its 

Tjlj- ' energy density) and there is sufficient reason to assume an even distribution of it over the space [for details see 

ref [6]]. Although DM energy density is expected to decrease at a faster rate than the density of DE throughout 
' the evolution, interestingly they have comparable magnitude today.. This surprising matching is known as 

the 'coincident problem'. To resolve this problem, use of tracker fields [7] and oscillating DE models [5| are 
normally employed. But recently, there arises a third possibility |9l [TOl [TTl [T^ by introducing DE and DM 
interact through an additional coupling term in the fluid equation. In the present work, we choose the third 
possibility as a solution of the problem. 

To have some inside about the unknown and mysterious nature of DE, many people have suggested that DE 
■ should be compatible with Holographic principle , namely "the number of relevant degrees of freedom of a system 

dominated by gravity must vary along with the area of the surface bounding the system" [13] . Such a DE model 
is known as Holographic DE(HDE) model. Further the energy density of any given region should be bound by 
that ascribed to a Schwarzschild black hole(BH) that fills the same volume jT4l[T5]. Mathematically, we write 

PD < MpL^'^, where po is the DE density, L is the size of the region(or infrared cut off) and Mp = (SttG) '^ 
is the reduced Planck mass. Usually, the DE density is written as 

PD = -JJ- (1) 

Here the dimensionless parameter 'c^ ' takes care of the uncertainties of the theory and for mathematical conve- 
nience the factor 3 has been introduced. In HDE paradigm [T31[Ill[Tni[I71[T51[Tn]one determines an appropriate 
quantity to serve as an IR cut off for the theory and imposes the constraint that the total vacuum energy in the 
corresponding maximum value must not be greater than the mass of a BH of the same size. By saturating the 
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inequality one identifies the acquired vacuum energy as HDE. Although the choice of the IR cut off has raised 
on discussion in the literature [151 EJ [201 ETJ [22j [23] , it has been shown, and it is generally accepted, that the 
radius of the event horizon of the universe(i?£;) the most suitable choice for the IR cut off where Re is defined 
as [H]: 

L = RE=ar-. (2) 
Jt 0, 

Now for the interacting DM and DE to resolve the coincidence problem (as mentioned above), the interaction 
term is chosen in the present work in the following two ways : (a) usually, the interaction term is chosen 
as AHpm + BHprt where pm and po are DM and DE densities and B are dimensionless constants. For 
convenience we shall choose A = B = 36^, i.e., interaction of the form Sb^Hp where p(= pm + Pd) is the total 
energy density, (b) a natural and physical variable interaction term is of the form "fPmPD with 7 a dimension full 
(mt ) constant. Note that this interaction term vanishes (as expected) if any one of the energy densities is zero 
while the interaction term grows with the increase of both the energy densities. Further, such an interaction 
term gives the best fit to observations [TT] for HDE models. In the present work we choose DM in the form of 
pure dust while the HDE, as perfect fluid with equation of state po = i-^dPd- 



2 Basic Equations 

We consider our universe to be homogeneous and isotropic flat FRW model and assume that it flUs with DM 
in the form of dust (having energy density pm) and HDE in the form of a perfect fluid having equation of state 
Pd = ^dPd where uje is variable. 

The Einstein field equations for spatially flat model are 

= p„ + PD (3) 

and 

2H=-p^-{1+ujd)pd (4) 

where for simplicity we choose SttG = 1 = c 

The conservation equations for the fluids are 

Pm + 'SHp„i Q (5) 

and 

Pd + 3H{1+ujd)pd = -Q (6) 

It is to be noted that for Q > 0, energy is transferred from DE to DM and opposite is the situation for Q < 0. 
As Q < would worsen the coincidence problem so we choose Q > throughout the work. Further, for 
validity of second law of thermodynamics and Lechatelier's principle [12 one must take Q > 0. Also it should 
be mentioned that baryonic matter is not included in the interaction due to the constraints imposed by local 
gravity measurements [24l [25] . 

Using the fleld equations Qand the acceleration of the universe is given by 

d = a(^H + H^^ =-^{p„, + {1+ojd)pd} (7) 

which shows that for the present accelerating phase it is necessary (but not sufficient) to have w/j < — i. 
Using the density parameters 

^m = {§^ , nn = ^ (8) 
the Einstein equation ([3]) can be written as 

+ nD = i (9) 
Introducing u = ^ as the ratio of the energy densities we have 

7; 1 
^m = ^— , no^^— (10) 
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3 Detailed Calculations for Q = 3b^Hp 



Using the conservation equations ^ and (jHl and the energy density of HDE from equation ([T]), we have the 
expression for the equation of state parameter as 



1 2^AI^ 



OJr) = 

3 



3c 



Qd ~3 



(11) 



which shows that there will be always acceleration. 
The evolution of the density parameter fi^i is given by 

fio = Hfll, (1 - Qd) 



1 
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and hence the ratio of the energy densities evolves as 



-u < 1 



cvTTd (1 — 



+ 36^ (1 + u) 



(12) 



(13) 



The Friedmann equation (|4]) and the conservation equation ^ can be converted(after a bit simplification) into 
an autonomous system as 

/ad 



Pd = '2pD 



H 



and 



H 



iH^ (1 - h^) 



PD 



ph 



3%/3c H 



(14) 



(15) 



The dynamical system has a line of critical points along the parabola po — 3c^H^ in the phase plane 
{pd, H) with the restriction 6^ = 1 — c^. Then the linearized matrix A has trace{A) = H (l — Ac^^ and 
determinant{A) — 0. So the phase paths form a family of parabolas [26j. The phase portrait for different 
choices of the parameter '5' are shown in figures 1(a) and 1(b). Note that along the line of critical points 

— and a;£) = — 1, u ~ Hence along the phase paths the ratio of the energy densities bears a 
constant value and the universe will be in the phantom era. 

Further, fixed points corresponding to u = is essentially a cubic equation which has at least one real root 
say, Uf. Then the parameter '6^' can be estimated by the fixed point as 
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1 + UfJ \3 3c^l + Uf 
Now, to analyse the stability of the fixed point we write 

2 



du du dt 
dx dt dx 



u 
H 



+ (1 + u) 



where x = In a. 

Then at the fixed point 



du' , 9 
— = 36^ 
du 



{u + 2) 



3{l + u)z 

1 2 - ^ 

1 + "/ c(l + u/)^ 



(16) 



(17) 



(18) 



U—U f 

< 0. 



Hence the fixed point is a stable one. 

Moreover the conservation equations ([5|) and ([6]) can be written as 



Pm (pm + pd) [b^PD - (l - 6^) Pri 



(19) 
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Phase Portrait for b=0.294 
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Fiff.l(a) 



Fig. 1(b) 



Fig. l(a)-l(b) represent the variation oipD—H. Though in Fig 1(b) the whole region is not a physically 
vaHd region but for better understanding about the system we have drawn the whole region. 



PD^- {pm+PD) [h^Pm + (l + WD + 6^) Pd] (20) 

From equation ([201) we see that pe, < 0, i.e., DE density decreases at least in the quintessence era. Also from the 
equation (119^ . if we assume |27| pu to be sufficiently large initially then matter density increases in the early 
phase and subsequently it decreases with pm = along the straight line pm — jz^Pd in the (p,„, po)— plane. 
Then in the phantom era (lud < —1), po may begin to increase and dominate over DM. 

Thus the present model of the universe shows a DE domination initially and subsequently the universe 
evolve with DM domination. Then there may be DE dominated phase at late time as predicted by observation. 
Hence this scenario is favourable for the present universe. 

Further, from the point of view of the coincidence problem we see that u is less than unity in the 

early phase of the universe and then it gradually increases, u ^ o(l) before Pm = or after p„i = or along 
the straight line pm — jz^Pd in the {pm, pd) -plane provided > or < or — \- Though the coincidence 

.2 

problem has partial solution around the straight line pm = jitipPd, but it does not give any explaination for 
u ^ o(l) in the present scenario. 



4 Calculation details for Q = "fPmPD 



Proceeding exactly as in the previous section the expression for equation of state parameter and the evolution 
of the density parameter are given by 



3 6c 3 



and 



tlD = HQd (1 - ^d) 



2x/tt; 



(21) 
(22) 
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Phase Portrait for 7=0.232,c=0.9 
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Fig.2(a) 



Fig.2(b) 



Fig.2(a)-2(b) represent the variation of ^Id—H. Here also the negative coordinate of H is not physicahy 
vahd. But for completeness of the system we have drawn the whole figure. 



So the evolution of the ratio of the energy densities is described as 



ii = 3Hu 



1 



7iJ 



3 3cVT 

Also the Friedmann equation (|4]) can be written as 



l + u 



H 



6 AC 



(23) 



(24) 



Hence equations <\22l and <\2A\ constitute an autonomous system in the phase plane {^D: H). The possible 
critical points are 
(?) = 0, is unrestricted, 
{ii) f^D = 0, ff = 0, 
(m) f2z3 = 1, H = Q , 
{iv)VLD — \, H is unspecified and 

(«) no = c\H^ ^. 

The first three critical points correspond to static model of the universe. In the first one the universe may have 
any amount of DE while in the second one the universe does not have any DE. The third and fourth critical 
points correspond to universe in phantom era filled only with the DE. For the fifth critical point, the ratio of 
the matter densities u — J and wd = — ^ <— 1. So the universe is again in the phantom region. A detailed 
investigation of this critical point will be done subsequently. 
The conservation equations ([5]) and ([B]) can be written as 



Pn 



IPD - ^/Mfi + Pd)p 



(25) 
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PD = -PD 



IPm + (1 + ^d) \/M~P m + Pd) 



(26) 



Apparently, we have similar situation as before, i.e., initially if we assume to have sufficient DE then p'o < 
and p'm > and subsequently p'm < 0. Note that p'm — along the curve ^"^ po^ — "iprr? — 'ipuPra' ■ But due 
to term in equation (j26p DE density begins to increase in phantom era and we should have DE dominated 
universe as expected in the present scenario. So both the energy densities are of comparable magnitudes (i.e., 
u ^ o(l)) twice during the evolution and give a possible explanation to the coincidence problem. 

To study the nature of the critical point ^c^, -i^^^ on the phase plane (S^b, H) we start with the linearized 
system : 



_ (l-c^) 



(2x 



y = 



(3-2c^)x 



3(l-c")y 
2c^ 



(27) 



where x ~ — c? , y = H ~ ^ . 
So the linearized matrix A has 



7 (l- c- 
tr{A) ^ L^-^ 
2 c^7 



> 



S 1 - c^) 
det{A) = --^^<Q 
2 c*7^ 

Thus for the linearized matrix A both the eigen values are real but of opposite sign and hence the critical 
point is of saddle type and unstable in nature. The phase portrait for different choices of 7 and are presented 
in figures 2(a) and (b). 

5 Discussion 

In the present work, we study the cosmological evolution of interactive DM and DE in the background of a 
homogeneous and isotropic FRW model of the universe. The DM is chosen in the form of dust while for DE we 
choose holographic DE in the form of perfect fluid having variable equation of state. The interaction between 
DM and DE is chosen either as a linear combination of the energy densities or in their product form. For both 
choices of the interaction term, the evolution equations can be suitably converted into an autonomous system 
and the critical points are analyzed both analytically and graphically. Finally, the present model of the universe 
shows partial solution of the coincidence problem. 
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